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ABSTRACT 


Suppose  one  observes  independent  Poisson  processes  with  unknown 

intensities  A.,  i  =  l,2,...,k,  and  that  apriori  it  is  believed  that  these 
c1  I 

intensities  satisfy  a  known  ordering.  For  preliminary  analysis,  it  might  be 
desirable  to  test  for  homogeneity  among  the  intensities  an  i,  of  course,  one 
would  want  a  test  that  utilizes  the  information  in  the  ordering.  hot  t. . 
denote  the  length  of  time  for  which  the  ith  process  was  observed. 

The  case  in  which  the  t.  are  equal  has  been  studied  in  the  literature.  We 
develop  the  conditional  likelihood  ratio  test  for  arbitrary  t_  .  This 

test  is  equivalent  to  the  unconditional  likelihood  ratio  test,  but  load:  to 
an  interesting  multinomial  testing  situation,  ie.  testing  for  homogeneity  o! 

p./t.  versus  a  trend  among  the  p./t.,  where  the  p.  are  the  cell  probabi 1 i t it 

11  li  l 

If  the  number  of  trials  in  the  multinomial  setting,  or  the  total  number  of 
occurrences  in  the  Poisson  processes,  is  large,  then  the  test  statistic  has 
approximate  chi-bar-squared  distribution  which  has  been  studied  in  the 
literature.  Results  of  a  Monte  Carlo  study  comparing  this  test  with  the 
maximin  test  developed  by  Y.  J.  Lee  are  discussed. 

Similar  results  are  also  obtained  for  testing  the  null  hy(x>thesis  that 
the  intensities  satisfy  the  prescribed  ordering. 
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1.  INTRODUCTION.  Barlow  et  al.  (1172)  discusses  the  problem  of  estimatin') 
a  finite  sequence  of  Poisson  intensities  which  are  assumed  to  be  noni nereas iriq . 
For  instance,  consider  a  system  which  is  observed  for  t  units  of  time  with  x 


1  . I 

failures,  is  then  modified  in  an  attempt  to  improve  its  performance,  is 


observed  for  t.  units  of  time  with  failures,  is  modified  aqain,  and  this 


is  repeated  until  it  observed  for  the  kth  time  for  t  units  of  time  with  X, 

k  k 


failures.  If  it  is  believed  that  the  modifications  will  not  harm  the  system's 
performance,  then  one  might  wish  to  estimate  the  vector  of  intensities, 


A  =  (A  , A  ,  .  .  .  ,  A,  )  ,  subject  to  A  >  \„  >  ...  >  A,  .  It  would  also  be  of 
12  k  1  —  2  —  —  k 


interest  to  test  for  homogeneity  among  the  intensities  with  the  alternative 


A  >  A  >  ...  A  and  A  >  A  ,  or  if  the  assumption  concerning  the  modification 

1  Z.  K  IK 


were  in  question,  one  could  test  A  A  ■  ...  >_  A  against  A.  -  >.  .  ,  tor  some  i 


i  +  i 


Suppose  X  , X.,...,X.  are  independent  Poisson  variables  with  means  ji  .  =  A . t . 

12k  i  l  J 


(0) 


let  be  a  partial  order  on  (1,2, ...,k),  let  A  be  a  fixed  vector,  let  a  be 


(0) 


an  unknown  scale  parameter  and  let  H  :  A  =  aA  ,  H  :  A.  •'  A.  whenever  i  j 

0  '  1  l  -  3 


and  (that  is,  A.  >  A.  for  some  i  <<  j).  The  hypothesis  H  ^  stipulates 


that  A  =  ( A  i  >  A  2 '  •  •  •  >  A  )  is  isotonic  (with  respect  to  <••)  and  we  suppose  that 


X  ^  is  isotonic.  We  consider  the  likelihood  ratio  test  (lrt)  for  11  versus 


rK 

H,  -  H_  and  H,  versus  H  conditional  on  )  .  ,X.  =  n.  While  it  will  ire  shown 

10  1  2  ''1  =  1  l 


chat  the  conditional  test  is  equivalent  to  the  unconditional  lrt,  it  does 
lead  to  an  interesting  multinomial  testing  situation.  W .>  also  know  that  for 
k  =  2  it  is  UMP  unbiased.  (See  p.  228  of  Ferguson  (1007).) 

Robertson  and  Wegman  (1078)  consider  order  restricted  tests  for  members 
of  the  exponential  family,  but  their  work  requires  that  the  sample  sizes  bo 
equal.  Their  results  can  be  applied  in  the  testinq  situation  considered  here 


only  if  the  t  are  all  equal.  Boswell  (1066)  considers  a  closely  related 


3 


problem  and  develops  the  conditional  lrt  for  testing  that  the  intensity  of  a 
nonhomogeneous  Poisson  process  is  constant  versus  it  is  nondecreasing. 

Lee  (1980)  developed  a  maximift  test  for  the  multinomial  setting 
arising  in  the  conditional  framework.  It  could  be  used  to  provide  a  test 
of  Hq  versus  -  H  .  The  results  of  a  Monte  Carlo  study  comparing  the 
lrt  and  the  maximin  test  are  given  in  Section  4.  It  is  found  that  for  small 
k  (k  =  3,4,5),  the  two  perform  similarly  with  the  power  of  the  maxinin  test 
larger  for  "regular"  alternatives  and  the  power  of  the  lrt  larger  for 
"nonregular"  alternatives.  Furthermore,  the  differences  in  power  are  not  too 
large.  However,  for  larger  k  (k  =  10)  the  differences  are  more  pronounced 
and  if  nonregular  alternatives  cannot  be  ruled  out  the  lrt  should  be  used. 

2.  ESTIMATION.  The  maximum  likelihood  estimate  (mle)  of  X  subject  to  11 
can  be  expressed  as  a  projection  onto  a  cone  of  isotonic  functions.  Wo 
introduce  some  notation.  With  <<  a  fixed  partial  order  on  {l,2,...,k}, 


let  denote  the  k-dimensional  reals,  let  C  =  {x  e  R^:  x  is  isotonic  with 
respect  to  <<},  let  w  =  (w  ,w  , ...,w  )  be  a  vector  of  positive  weights,  for 

X  x  K 

vk 

x,y  t.  R  let  (x,y)  denote  the  inner  product  2.-_iw-xY'  anci  f°r  V  ‘  **'* 

K  W  X.  —  1  XXX.  K 

Ew(y|c)  denote  the  projection  of  y  onto  C,  that  is  E^tylc)  minimizes 

yk  ,w. (y.-x.)2  for  x  e  C. 

Li=l  li  i 

Theorems  1.4  and  1.5  of  Barlow  et  al.  (1972)  state  that 

(1)  1*  . w . E  (y | C)  .  =  w.y . 

^1=1  l  w  '  1  4'1  =  1  11 


(2)  w  (y.-E  (y}c).)E  (y}c) .  =  0 

*1  =  1  i  l  w  '  1  w  1 

The  mle  of  X  subject  to  H  is  X  =  e  (X/t|c)  where  t  -  (t  ,t  ,...,t  ), 

x  t  1  x  K 

X  =  ^Xl'X2’ '  '  ‘ 'Xk*  and  f°r  X,y  X//y  =  (x/yi'X2/y2 . W 


xy  =  (x,y, ,X„y_, • • . .x  y  ) .  (See  Barlow  et  al.  (1972)  page  44.) 
1  1  2  2  k  k 


4 


Conditioning  on  Y  =  £^_^X^  =  n»  the  density  of  X  is  that  of  a  multi¬ 
nomial  with  parameters  n  and 

(3)  p.  =  X.t./Ik  ,X.t.  for  i  =  1,2, - ,k. 

l  li  ^3=1  3  3 

So  H  :  A  =  aX ^  is  equivalent  to  H':  p.  =  pf^  =  x|°^t  /) k  ^ I  ,  I  i  k , 

0  0  i  *i  i  i  '  j- 1  3  j 

H^:  X  is  isotonic  is  equivalent  to  11^:  p/t  is  isotonic,  and 
H^s  X  is  not  isotonic  is  equivalent  to  :  p/t  is  not  isotonic. 

The  mle  of  p  subject  to  is  also  of  interest. 

Theorem  1.  The  mle  of  p  subject  to  is  given  by  p  =  tE^ (X/t | C) /n.  These 
p.  satisfy  p.  >_  0  for  1  <_  i  k,  £k_^p.  =  1,  and  p  -+  p  almost  surely  provided 
p/t  is  isotonic. 

k  Xi 

Proof.  The  mle  of  p  under  H’  maximizes  n.  , (p./t.)  subject  to  p/t  is 
-  1  1=1  11 

isotonic,  p.  >  0  and  Tk  ,p.  =  1.  Applying  the  result  on  p.  46  of  Barlow 

l  —  ‘’1=1  l 

et  al.  (1972),  we  see  that  p  =  tE^ (X/t| C) /n.  Since  the  projection  E  ( * | C) 
is  continuous  for  fixed  C  and  w,  p  =  tEfc (X/ (nt) j C)  ->  tEt(p/t|c)  almost  surely. 
If  p/t  is  isotonic  the  right  hand  side  (rhs)  is  p  and  the  proof  is  completed. 

If  <<  is  a  total  order,  then  the  pool-adjacent-violators  algorithm  can 
be  used  to  compute  the  projection  in  the  formula  for  p  and  the  lower  sets 
algorithm  can  be  used  for  an  arbitrary  partial  order.  (See  Chapter  2  of 
Barlow  et  al.  (1972).) 

3.  TESTS  OF  HYPOTHESES.  As  was  mentioned  in  the  Introduction,  we  consider 

the  conditional  lrts  of  versus  H  -H^  and  °f  versus  H^-  However,  these 

lead  to  the  lrts  of  H'  versus  H'  -  H'  and  of  11'  versus  H' .  Chacko  (1966) 

0  10  1  2 

developed  the  lrt  and  an  asymptotically  equivalent  modified  test  for  IP 
versus  in  the  totally  ordered  case  with  p^  and  t  constant  vectors . 

Robertson  (1978)  developed  both  lrts  for  partial  orders  with  t  a  constant 


vector. 


We  se<  (hat 


Denoting  the  conditional  likelihood  ratios  by  A^  and 

(0) 

<4>  TM  =  -2  ln  Ki ‘  2tixillnEt't|c,i  -  ln  - 2n  lE " 

1 

and 

(5)  TJ  -  -2  1»  X-  -  2l^1X.(l„  a.  -  J.  Et<||c)  .  I. 

1 

Hence,  the  test  statistic  for  versus  n  -  Hq  is  obtained  by  replacing 

pf^  with  A(0)t./Vk  ,  A^^t.  in  (4)  and  for  testing  H,  versus  H  ,  T  =  T'  . 
l  li  1 j  =  l  3  j  1  2  J  2  1 2 

Remark  2 .  The  conditional  lrts,  T  and  T  ,  are  also  the  unconditional  Irts. 

Proof.  Under  H  ,  the  mle  of  a  is  n/Y^  ,  A^^t  .  Using  (1)  and  straight-forward 
—  0  ']=1  3  3 

algebra  and  denoting  the  likelihood  ratio  for  testing  11^  versus  H^  -  II 

(H  versus  H  )  by  A  (A  )  ,  one  can  show  that  -2  ln  A  =  T  and  -2  In  A  -  T  . 

J.  2  UI  x4  U 1  u  I  i 

Next  the  large  sample  distributions  of  these  test  statistics  are 
determined.  The  derivations  are  like  those  given  in  Robertson  (1U78)  and 
in  fact  we  will  make  use  of  several  lemmas  proved  there. 

Lemma  3.  In  the  multinomial  setting,  let  "  be  independent  normal  variables 

_  ^ 

with  mean  0  and  variance  1/p.  and  let  Z  =  )  .  , p.Z..  As  n  4  ®, 

■  V  _  _*  J’1  1  1 

t'rTip-p)  >  (p  (Z  -Z)  ,  p  (Z  -Z) ,  ....  p  (Z  -Z)  )  .  bet  q  denote  a  fixed  vector 

IX  2  2  k  k 

of  probabilities  with  q/t  isotonic,  let  D  denote  the  closed,  convex  cone 


{x  e  R x.  <  x.  if  i  <<  j  and  q./t.  =  q./t.},  let  r,  <■  r_ 

k  l  -  3  1133  12 


r,  denott 
h 


the  distinct  values  among  q^/t^,  i  =  l,2,...,k,  and  let  M(i)  =  ij:  cj ^ / 1. ^  -  r^ 

for  i  =  l,2,...,h.  Robertson  observed  that  E  (x|d)  can  be  computed  by 

w 

independently  computing  its  values  for  subscripts  in  M(i)  for  i  =  l,2,...,h. 


This  implies  that  if  y  t  has  positive  entries  and  is  constant  on  each  M(i) ,  then 

(6)  E  (x  |  D)  =  E  (x|d)  . 

w  wy 

(Tliis  can  be  easily  seen  by  considering  the  lower-sets  algorithm.) 

Lemma  4.  If  x,y  1  R,  with  y  constant  on  each  M(i),  then  E  (x-y  D)  - 
-  k  w 

E  ( x  j  D)  -  y.  If  in  addition,  y  has  nonnegative  entries,  then  E  (xy|D)  -  yE  (x|l>). 
w  w  w 
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Lemma  5.  If  x  c  R  with  max.  .....  x.  <  min.  .. . .  x.  for  i  = 

-  k  ] cM ( 1)  1  3cM(i+1)  j 

then  E  (xlc)  =  E  (x|d) . 
w  w  1 

The  cone  D  is  determined  by  <<  and  q.  For  a  given  D,  define  P  (H,k)  to 

w 

be  the  probability  of  exactly  i  distinct  values  in  E  (v|d)  where 

V  =  (V, ,V.  . . . ,V,  ),  with  the  V.  independent  normal  variables  with  mean  0  and 
12k  i 

variance  1/w. .  If  <<  is  the  usual  total  order  on  fl,2,...,k)  and  q  /t  is 
r  i  i 

constant,  then  D  =  {x:  x.  <  x  „  <  ...  <  x  )  and  the  P  (t.,k)  are  discussed  in 

1  —  2  —  —  k  w 

detail  in  Barlow  et  al.  (1972) .  Approximations  for  nonconstant  w  are 

discussed  in  Siskind  (1976)  and  Robertson  and  Wright  (1983).  If  is  the 

same  total  order,  but  q./t.  is  not  constant,  then  the  results  mentioned  above 

can  be  used  to  determine  these  probabilities  for  each  M(i)  and  the  h-fold 

convolution  gives  the  desired  P  U,k).  (cf.  p.  148  of  Barlow  et  al .  (1972)). 

w 

The  latter  reference  also  discusses  the  P  (£,k)  for  some  partial  orders. 

w 

2  2 
Let  Xj£  denote  a  chi-squared  variable  with  i  degrees  of  freedom  (xQ  0)  . 

Theorem  6.  If  in  the  Poisson  setting,  the  vector  of  intensities  is  a  X  ^ 

and  q  =  p^  (cf.  (3)),  or  in  the  multinomial  setting,  if  the  probability 

vector  is  q  =  p*°^  and  if  D  is  determined  by  q,  then  for  t  ^  0, 

(7)  lim  _1?[t  >tjY=nl  =  lim  p[t'  >t]  =  V*  ,P  (S.,k)P|v2  >tl  . 

nr>®  1  01—  1  J  n-^  1  01—  1  Ll=l  q  1  *5,-1—  J 

If  in  the  Poisson  setting,  the  vector  of  intensities  is  of  the  form  av  with 

v  isotonic  and  q  is  determined  by  (3)  with  X  ^  replaced  by  v,  or  if  in  the 

multinomial  setting  the  vector  of  probabilities,  q,  is  such  that  q/t  is 

isotonic,  and  if  D  is  determined  by  q,  then  for  t  >  0, 


(8) 

lim  P [T  >t 1 

n-*»  1  12-  1 

Y=nJ  =  limn^,PfTi2lt]  =  ^=1pq(^k)P|X 

(9) 

lim  P[T,  >t[ 

n-**>  12—  1 

Y=n I  £  limn wPe[Tl2>t|Y=nf,  and 

(10) 

lim  P[T*  >t] 

n-*°°  '  i2-  1 

<  lim  P,[T'  >t], 

-  n-w  «Sl  12-  J 

where  in  the  Poisson  setting  P  |)  denotes  the  probability  under  X  -  (1, !,...,!) 


7 


and  in  the  multinomial  setting,  P  []  denotes  the  probability  under 


5.  =  t./y .  . t . , 

l  i  '']=1  3 


1,2, . . . , k . 


Comments .  In  the  Poisson  setting,  if  one  wishes  to  test  H  versus 

(0) 

,  then  q  is  set  equal  to  p  determined  by  (3).  The  vector  q 
determines  D  which  in  turn  determines  the  P  (£,k)  and  so  large  sample 
p-values  can  be  calculated  for  the  conditional  test  from  (7) .  In  testing 
versus  H^,  (9)  indicates  that  the  asymptotically  least  favorable 

configuration  in  H  is  X  =  (1,1,. ..,1)  and  so  with  q.  =  t./)k  t., 

1  li  '3=1  3 

D  =  C  and  approximate  p-values  can  be  computed  from  (8) . 

In  the  multinomial  setting,  in  testing  versus  -  H^,  set 

(0) 

q  =  p  and  use  (7)  to  determine  large  sample  p-values.  In  testing 
versus  ii^,  the  asymptotically  least  favorable  configuration  in  H  is 
q^  =  t./^t.  in  which  case  D  =  C  and  approximate  p-values  can  be 
computed  from  (8) . 

In  either  case,  one  might  not  want  to  use  the  asymptotically  least 

favorable  configuration  when  testing  versus  H ^  or  versus  so  one 

could  use  the  restricted  mles  of  X  or  p  rather  than  e  or  6. 

Proof.  T_,  is  defined  to  be  T'.  with  p^  replaced  bv  X^t  X^t 

-  01  01  1  i  1  ' '  j  =  1  j  3 

and  T  =  T'  .  Furthermore,  conditional  on  Y  =  n,  (X  .X  ,...,X, )  is 
12  12  12k 

multinomial  with  parameters  n  and  p.  =  X  t  /Yk  X  t,  and  H  li! ,  i  =  0,1, 

1  11  3=1  33  11 

So  we  need  only  consider  the  multinomial  situation.  We  first  consider 


the  distribution  of  T^  under  H^.  Setting  p^  =  xVn,  expressing  as 

(0) 


2nyi=1Pifln  E't(p/tjc)i  -  lnfp^'/tj}. 


and  expanding  In  E(p/t|c)  .  and  ln(pf°Vt.)  about  p./t.,  we  write  T'  as 
t  1  1  11  11  01 


2nIi=lti(Et(p/tlc)i  ~  pi0)/ti} 


+  nX^_1Pi ^ ( (pi“ pi°^ ) /fci ) ~  (Et<p/tlc)i  “  Ib/tj) » 


2  ,,,2 


2,2, 


(11) 


8 


where  a.(B.)  is  between  p./t.  and  E  (p/t| C) . (p!° Vt. ) .  Under  H',  both 
11  11  t'l  ii  0 

E^_  (p/t |  CK  and  pVK  converge  almost  surely  to  p^°^/t  .  Recall,  (1) 

implies  that  the  first  term  in  (11)  vanishes. 

Since  p/t  is  consistent  for  p^/t  under  11^,  there  is  for  almost  all 

uj  in  the  underlying  probability  space  an  N,  possibly  depending  on  u>,  for 

which  p/t  satisfies  the  hypothesis  of  Lemma  5  for  n  >_  N.  Hence,  for  such 

a)  and  n.  Lemmas  4  and  5  can  be  applied  to  the  second  term  in  (11)  to  obtain 

.(OK  .  „  ,  r,;  (OK  ,  . .  (0) 


</n(p.-p.  )  t.E  (t/n(p-pl  )  / 1 1 D)  .  -  4i(p.-p.  ,  , 

- 1  — — )  -  (-^ - - - - - — )  } , 


) 


t.B. 
i  i 


t.a. 
l  l 


which  converges  in  distribution  to 


Ic  fc2  P-  ,  (0)  p. 

h-i  1  <-t“<Vz> >  -  IEt'V(z-z|lDli  -  -r-,zi- 


(0) 


(0) 


(0) 


-n2, 

^ )  )  ) . 


Pi 


Applying  Lemma  4  and  (6) ,  this  can  be  written  as 


(12) 


Ik  lP(0)(z.-z)2  -  lk  iP^(Z.  - 

1=1  1  1  u 1=1  1  x 

■ «  (o,(ziD,i-z)2  *  2tipri,zrE  (o)<ziD,i,<E 


E  (0) <2lD)±> 

p 


(0)  (r‘lD)i-Z) 


The  second  term  in  the  rhs  of  (12)  can  be  shown  to  be  zero  using  (1)  and 
(2),  and  Theorem  3.1  of  Barlow  et  al.  (1972)  shows  that  the  first  term  on 
the  rhs  has  the  desired  distribution. 

Next  we  consider  the  distribution  of  T^  with  q  a  fixed  probability 
vector  for  which  q/t  is  isotonic.  Writing  T^2  as  2n£^  p^{ln(pVtJ -In  E^fp/tlo 
expanding  In  Et(p/t|cK  about  pVK  and  applying  (1),  Tj^  can  be  written  as 


tlPi(^(Et<P/tlC)i  "  Pi/ti))2/Yi 

where  y^  is  between  E^(p/t|cK  and  p/t  and  hence  converges  almost  surely 
to  q^/t  .  As  in  the  proof  of  the  first  part  of  the  theorem,  for  almost  all 
w  and  n  sufficiently  lar  'e  /n (F^ (p/t | C) -q/t)  =  E  ( /n (p-q) /t | D) . 


'3 


Hence,  T|2  converges  in  distribution  to 


t2 


q. 


(13)  £  —  (E  (^<Z-Z)  |D)  .  -  ~(Z.~Z))Z  =  a  (E  (z|d).  -  Z.)2. 

q  t  t  1  i  t  i  ^i-1  i  t  i  i 

i  i 

Applying  (6)  this  becomes 

lk  a. (E  (z|d) .  -  z. ) 2 

1=1  i  q  1  i  l 

and  this  has  the  desired  distribution  (cf.  Theorem  2.5  of  Robertson  and 
Wegman  (1978) . ) . 


We  establish  (10)  to  conclude  the  proof.  The  variables  U.  =  fu./t.  . 

1111 

are  independent  normal  variables  with  means  zero  and  var(U.)  =  1/t,.  Using 

i  i 

(6),  the  rhs  of  (13)  can  be  written  as  Vk  ,t.(E  (u|d).  -  U.)2,  which  is  the 

^1=1  i  t  l  i 

distance  from  E  (u|d)  to  U.  Since  CCD,  by  the  definition  of  projection 

this  is  maximized  for  D  =  C,  which  occurs  if  q/t  is  constant,  ie.  q.  =  t./V^  .  t . . 

li  “i=l  j 

In  the  multinomial  setting,  Lee  (1980)  developed  a  maximin  test  for 

Pl/tl  =  P2/fc2  =  =  versus  Pi+1/fci+1  XdpVt^  i  =  l,2,...,k-l,  with 

rk 

d  >  1.  The  test  statistic  is  ^ix^,  which  has  an  approximate  normal 

distribution  and  under  this  null  hypothesis  its  mean  and  variance  are 


n n  ,it./r  .t.  and  n{[k  ,t.  -  ([*  .it./^  .t.)2}, 

‘■1=1  l  ^i=l  l  ^i=l  i  ^1=1  i  Li=l  i  ‘-i=l  i 

respectively.  The  tests  S  and  T^  are  compared  in  the  next  section. 

4.  COMPARISON  OF  THE  LRT  AND  MAXIMIN  TESTS.  The  maximin  test  is  a  contrast 
test  and  Section  4.2  of  Barlow  et  al.  (1972)  contains  a  discussion  of  the  use 
of  the  likelihood  ratio  and  contrast  approaches  in  testing  for  trends  among 
normal  means.  They  concluded  that,  while  the  contrast  test  is  typically 
much  easier  to  use,  the  lrt  provides  the  most  satisfactory  general  way  of 
incorporating  prior  information  about  ordering.  In  the  case  of  a  total  order 
and  small  k,  the  contrast  statistic  provides  a  suitable  alternative.  If 
additional  information  is  available  about  the  spacing  of  the  parameters,  then 
a  contrast  test  based  on  this  additional  information  may  be  preferred. 


To  give  some  idea  of  the  differences  in  the  power  for  the  two  to;, to  ,1 

Monte  Carlo  study  was  conducted.  In  particular,  if  one  were  testing 
rk 

Pi  =  ti/Zj=1tj'  versus  Pj/t^/t^. •  wiLh  i'1/tJ ' v/\' 

then  T’  ,  with  pf°^  replaced  by  t./)k  t.,  and  S„,  could  be  used.  because 

01  x  i  '-3=1  3  01 

the  distribution  of  under  H  ^  is  quite  complex,  Monte  Carlo  exper  imer.ts 

were  conducted.  With  k  =  3,  4,  5,  10,  n=25,  80,  nominal  levels  of  .1,  .05, 

and  various  choices  of  t  and  p,  the  powers  of  T'  and  5;  ,  were  approx i  mated 

01  0l 

based  on  5000  repetitions.  Some  of  these  values  are  given  in  Tables  1  and  2. 

To  assess  the  accuracy  of  the  approximations  for  the  distributions  under  the 
null  hypothesis,  the  estimates  of  the  power  under  the  null  hypothec, i are 
included.  For  n=80  and  the  cases  presented  in  these  tables  the  largest; 
discrepancy  in  the  a  level  for  T^S^)  and  a  nominal  level  of  .1  is  .Ole  ( .  uo  ’ ) 
and  for  a  =  .05  it  is  .005  (.007).  For  n  =  25  the  maximum  discrepancies  were 
larger  but  both  approximations  seem  to  be  useful  for  k  in  this  range. 

However,  for  k  =  10  and  n  =  25  the  approximation  for  the  distribution  of 
seems  to  give  a  test  that  is  quite  conservative.  Its  estimated  a  level  is  .005 
(.031)  when  the  nominal  level  was  .1  (.05).  The  approximation  for  S  seemed 
quite  adequate  even  with  n  =  25  and  k  =  10.  Its  estimated  a  levels  are 
.100  (.046),  respectively. 

It  is  clear  from  Tables  1  and  2  that  neither  test  is  uniformly  better 
than  the  other.  In  fact,  when  the  p/t  increase  regularly,  such  as  in  the 
cases  p/t  =  (.25, .30,. 45) ,  p/t  =  ( . 20, . 25, . 35) ,  p/t  =  ( . 15 , . 20 , . 30 , . 35) , 
p/t  =  ( . 10, . 15, . 20, . 25, . 30) ,  etc.,  then  the  maximin  test  outperforms  T^  , 
but  for  irregular  increases  in  p/t,  such  as  in  the  cases  p/t  -  ( . 25 , . 25 , . 50) , 

P/t  =  (.2, .3, .3),  p/t  =  (.052, .052, .120) ,  p/t  =  ( . 04 , . 04 , . 04 , 07  i )  ,  p/I 
( . 03 , . 03 , . 07 , . 07) ,  p/t  =  ( . 15 , . 2 , . 2 , . 2 , . 25) ,  etc.,  T^  has  qreater  power 
than  S^.  For  k  =  3,4,5  the  differences  in  power  are  not  too  large  and  the 


magnitudes  are  similar  in  both  directions.  So  for  small  k  one  could  us<- 
S  ^  if  the  alternative  were  believed  to  be  "regular"  in  the  sense  described 
above  or  T^  could  be  used  if  it  is  desirable  to  protect  against  nonregular 
alternatives . 

It  is  interesting  to  note  that  for  k  =  3  and  4  the  above  conclusions 
held  whether  the  vector  t  was  constant  or  not.  (Several  other  choices  of 
p  and  t,  not  given  in  Tables  1  and  2,  were  considered  and  these  conclusions 
were  substantiated  in  those  cases,  also.)  For  this  reason  only  const  a:.*,  t  ' 
were  considered  for  k  -  5  and  10.  Recall  that  for  k  5,  the  1-  ,ki 
are  intractable  for  nonconstant  t. 

Power  comparisons  were  made  for  k  =  10,  but  with  n  25  they  were  not 
very  meaningful  because  of  the  conservative  nature  of  the  approximation  to 
the  null  distribution  of  .  For  n  =  80  both  approximations  were  very 
reasonable  and  so  power  comparisons  could  be  made  in  that  case.  Linear 
alternatives  were  considered.  For  =  i/55  both  tests  have  powers  that 
are  essentially  one.  So  p.  =  i/110  +  .05  was  considered.  The  teste,  with 
nominal  level  .1  had  powers  .822  and  .861,  and  the  tests  with  nominal  1 ( vc 1 
.05  had  powers  .712  and  .759.  Of  course,  the  maximm  test  performed  U  1  ter 
for  such  an  alternative.  The  alternative  =  .09,  1  i  /  9,  and  p^  .1 
was  considered.  The  estimated  powers  for  the  tests  with  a  =  .1  are  .700 
and  .556  and  for  the  tests  with  a  =  .05  they  are  .578  and  .409.  Finally  th 
alternative  =  .07,  p^  =  ...  =  p^  =  .1,  and  =  .13  was  considered,  and 

the  approximate  powers  for  the  a  =  .1  tests  are  .377  and  .332.  For  the  <  - 
tests  they  are  .252  and  .199.  Of  course,  the  last  two  alteri itivos  are 
nonregular  and  the  lrt  has  the  larger  power.  In  these  cases,  (he  i  ii>-r«  •.  t  .■  ;< 
in  power  may  be  as  large  as  40  percent  and  so  the  lrt  should  definitely 
be  considered  to  guard  against  nonregular  alternatives  for  larger  k. 
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TABLE  1.  ESTIMATED  POWERS  OF  THE  MAXIMIN  AND  LRTS ,  k=3 


k=3 

n=25  n=80  n=25  n=80 

Nominal  Level 


.10 

.05 

.10 

.05 

.10 

.05 

.  10 

.05 

t=(l, 

.1.1),  p/t= 

d/3,1/3, 

1/3) 

t= ( 2 , 1 , 

,  1) ,  p/t= 

(.25, . 

.25, 

.25) 

T' 

01 

.084 

.059 

.108 

.047 

.101 

.047 

.085 

.052 

*01 

.088 

.056 

.095 

.043 

.081 

.049 

.  100 

.044 

t=(l, 

.1.1),  p/t= 

(.25, .30, 

-45) 

t= (2,1, 

.  1) ,  p/t= l 

(-20, . 

,25, 

•  35) 

Toi 

.389 

.227 

.718 

.565 

.459 

.  315 

.771 

.664 

S01 

.451 

.350 

.812 

.689 

.423 

.  337 

.802 

.677 

t=(l,l,l),  p/t=(.25, .25, .50)  t= (2,1,1) ,  p/t= ( . 2 , . 3 , . 3 ) 


T* 

01 

.559 

445 

.937 

.871 

.346  .220 

.624 

.509 

.572 

475 

.917 

.843 

.304  .230 

.633 

.473 

t=(2, 3, 5) 

,  P/t= 

(.10, .10, .10) 

t= ( 2, 3 , 5) ,  p/t= 

(.05,. 

10, . 12) 

T* 

01 

.126 

046 

.093 

.054 

.571  .332 

.89) 

.803 

.095 

056 

.106 

.048 

.504  .386 

.889 

.778 

t=  ( 8 , 2 , 8 ) 

,  p/t= 

(1/18,1/18,1/18) 

t= (8, 2 , 8) ,  p/t= 

(.04,  . 

06, .07) 

T0l 

.112 

051 

.084 

.052 

.527  .351 

.832 

.754 

S01 

.080 

053 

.100 

.051 

.452  .368 

.855 

.763 

t=(4,6,4) 

,  P/t= 

(1/14,1/14,1/14) 

t= (4 , 6 , 4 ) ,  p/t= 

(.052, 

.052, .120) 

T' 

01 

.112 

067 

.101 

.048 

.758  .609 

.981 

.947 

.117 

040 

.100 

.046 

.723  .535 

.965 

.921 
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TABLE  2.  ESTIMATED  POWERS  OF  THE  MAXIMIN  AND  LRTS,  k=4,5 


k=4 


n=25  n=80  n=25  n=80 

Nominal  Level 


.10  .05 

.10 

.05 

.10 

.05 

.10 

.05 

t=(l, 1,1,1) , 

p/t=( .25, 

.25, .25, .25) 

t= (4,6, 

4,6)  , 

p/t= ( . 05 , 

.05, .05 

,.05) 

T01 

.089  .055 

.097 

.051 

.099 

.054 

.095 

.05  5 

S01 

.105  .050 

.101 

.050 

.086 

.042 

.  103 

.049 

t=(l, 1,1,1) , 

p/t= ( . 15 , 

.20, .30, .35) 

t=  ( 4 , 6 , 

4,6)  , 

p/ t= ( . 034 

,  .042,  . 

0495,  . 

Toi 

.578  .446 

.930 

.860 

.487 

.338 

.826 

.729 

S01 

.634  .482 

.945 

.883 

.480 

.342 

.860 

.757 

t=(l, 1,1,1)  , 

p/t=(.10. 

.25, .25, .40) 

t=(4,6, 

4,6)  , 

p/t=( .04, 

.04, .04 

,  .07  3) 

T01 

.752  .665 

.995 

.985 

.501 

.345 

.844 

.766 

S01 

.804  .668 

.993 

.981 

.464 

.333 

.834 

.727 

t-(l, 1,1,1) , 

p/ 1= ( . 2  3 , 

.23, .23, .31) 

t= (4 ,6, 

4,6)  , 

p/t=(.03, 

.03, .07 

,  .07) 

T01 

.208  .129 

.376 

.256 

.730 

.562 

.983 

.955 

S01 

.240  .143 

.388 

.245 

.684 

.544 

.981 

.946 

k 

=5  and  t=(l. 

1,1,1, 1) 

p/t=(.2, 

.2, .2, .2, . 

2) 

p/t=(. 

10,  .1 

5, .20, .25, 

.  30) 

Toi 

.094  .050 

.099 

.049 

.649 

.489 

.  9b6 

.930 

soi 

.086  .051 

.091 

.050 

.679 

.567 

.974 

.946 

P/t= ( . 15,  . 

15, .15, .15 

,.40) 

P/t=  ( 

.15,. 

2, .2, .2, .25) 

T01 

.711  .613 

.985 

.969 

.  264 

.172 

.496 

.  3  50 

b01 

.654  .555 

.955 

.924 

.252 

.176 

.488 

.  361 

069) 


